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Abstract. We reanalyse Cosmic Microwave Background data from experiments probing 
ly-^ , both large and small scales. We assume that measured anisotropies are due not only to 

primary fluctuations but also, especially at small scales, to secondary effects (namely the 
\Q , Sunyaev-Zel'dovich effect) and possible point source contaminations. We first consider pri- 

o 

mary and secondary anisotropies only. For the first time in such analyses, the cosmological 

rG ' ' ' 

dependence of secondary fluctuations is fully taken into account. We show in that case that 

q . a higher value of the normalisation erg is preferred, as found by previous studies, but also 

■ higher values of the optical depth t and power spectrum index n s are needed. In the sec- 

CZJ . 

ond part of our analysis, we further include possible contaminations from unresolved and 

unremoved point sources. Under these considerations, we discuss the effects on the cosmo- 
■ *~ i ■ 

' logical parameters. We further obtain the best combination of relative contributions of the 

\—{ ' 

three kinds of sources to the measured microwave power on small scales at each frequency. 
Our method allows us to simultaneously obtain cosmological parameters and explain the 
so-called small scale power excess in a consistent way. 
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1. Introduction 

Measurements of the Cosmic Microwave Background (CMB) angular power spectrum are now 
becoming invaluable observables for cosmology. The detailed shape of this spectrum allows 
one to determine cosmological parameters with high precision. The CMB anisotropy spec 



trum has been recently measured from large to intermediate scales with WMAP (Benn ett et al 



Send offprint requests to: Marian Douspis 



Marian Douspis et al.: Small scale contributions to CMB: A coherent analysis 



2003), Archeops (Benoitetal 



2003) and BOOMERanG (Jones etal 



small scales, differen t experi ments, ACBAR I Kuo et alJ 120 04). CBI (Pearson etal 



Readhead et al 



2004a), VSA (Reb olo et al 



34) 



2005) expe riments. At 



2004) and BIMA ( Dawson et al 



2003 



2002), have de- 



tected signal in excess of that expected from purely primordial CMB anisotropies. This excess 
is not fully un derstood yet but several explanations have been proposed. Some relate to the early 
Universe (e.g. Cg^r^y_^Mgl chiorril2002 ; Griffiths et al 12003 ), some are associated with reionis- 
J2003I) " 



ing sources (e.g. 



Oh et al. 



and others are related to the astrophysical processes contributing 



to the signal at the CMB frequencies (e.g. 



Bond et al. 2005 



Toffolatti et al. 2005). As a matter of 



fact, when observing at microwave wavelengths one expects to detect not only primordial CMB 
fluctuations but also secondary anisotropies as well as galactic and extra-galactic foregrounds. 
Astrophysical contributions are, as much as possible, removed during data analysis to produce 
clean estimates of the CMB angular power spectrum. Nevertheless, residuals may remain and 
contribute to the measured power. 

In order to minimise galactic contamination, most of the small sky coverage experiments 
(probing small scales) observe at high galactic latitudes in regions where the sky is as much as 
possible free of galactic emissions. Moreover, the galactic signal can be removed through multi- 
frequency observations, or even by purely masking the contaminated regions. Since most of the 
residual signal is expected on large scales, galactic emissions only very moderately affect small 
scale CMB observations. 

At small angular scales, additional contribution to the primary CMB signal arises from the 
interaction of CMB photons with matter along the ir propagation path. These are the so-called 
secondary anisotropies (e.g. 



!Iu & Dod jlsonl200 2. and references therein). The lader are dom- 
inated by the thermal Sunyaev-Zel'dovich (SZ) effect 1 Sunvaev & Zerdovichlll972 , 



1980), i.e 



CMB photon inverse Compton scattering off free electrons in the hot intra-cluster medium. When 
clusters are known and resolved by an instrument, they can easily be cleaned out of CMB maps. 
However, not all clusters are known and we expect a contribution to the microwave signal from 
unresolved/unknown clusters. The SZ effect being frequency dependent (contrary to primary 
CMB fluctuations) and its spectral signature being known, one can use this information to detect 
and remove SZ signal from unknown clusters in a map if multi-frequency observations are con- 
ducted. On the contrary, with single frequency experiments, we are consequently left with an SZ 
contamination at small scales. 

Contamination by extra-galactic sources has long been studied for CMB and SZ ob- 



servations. The effects of radio sources 


were investigated for low frequency experiments 


(e.g. 


Ledlow & Owen 


1996 




Coorav et al. 






Toffolatti et al. 


1998 




Sokasian et al. 


2001: 


de Zotti et al. 


2QQ5 


Gonzalez-Nuevo et al. 


122j 

2QQ5 


). At high frequencies dusty galaxies emitting 



Blain 1998 



White & Maiumdar 



in the infra-red are the dominant sources of contamination (eg 
2004 ; La gache et alJ2005 ). Contamination by point sources can be monitored by interferometric 



arrays, and multi-frequency observations should help us in reducing the additional signal through 
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component separation. Once again, however, a residual contamination is likely to remain, espe- 
cially in view of the uncertainties in the modelling of the sources themselves. 

The excess power reported by small scale CMB experiments may thus be due to several 
astrophysical signals expected to contribute at those scales. In order to correctly estimate the 
relative and possible contribution of the secondary effects or extra-galactic signal on top of the 
primary CMB anisotropy, one needs a coherent analysis. Cosmological dependences are expected 
for CMB and SZ effect, and SZ and point sources are contributing o n rather similar scales. Thus, 
unlike previous studies (e.g. 



Goldstein et al. 2003 



Bond et al. 2005), we do not decorrelate small 



scales from those on large scales in our analysis of the CMB data. On the contrary, we look for 
the best model that fits data simultaneously on large and small scales. Hence, assuming that 
primary, secondary and point sources signals are all contributing to the model, we present a 
coherent analysis (i.e. including the full cosmological dependence of secondary effects) of the 
CMB angular power spectrum as measured from large (£ = 2) to small (i = 10000) scales. 

In the next two sections we describe the CMB data we used and we present our modelling of 
the power spectra of the different contributions that enter the analysis. We then present our results 
in terms of cosmological parameters inferred from an analysis accounting coherently for primary, 
secondary and extra-galactic contributions. We show in each case the relative contributions of all 
considered signals, and discuss our results. 



2. Observations 



CMB temperature anisotropics have been now measured on a wide range of scales with high sig- 
nificance. WMAP has yielded the best measurements from large to intermediate angular scales, 
and together with Archeops and BOOMERanG, offers an estimate of a series of features (peaks 
and troughs) consistent with those expected from acoustic oscillations in the primordial fluid. 
ACBAR, CBI, VSA and BIMA observed smaller scales, those affected by exponential damping 
of power of primary fluctuations due to photon diffusion and to the finite thickness of the sur- 
face of last scatterin g. The C MB polarisation (TE, EE, B B) has also been measured by several 
exper iments (£DASI, 



2005), (CBI, 



Leitch 



et al 



Readhead et al 



2005), (BOOMERanG, 



2004b), (WMAP, 



Piacentini et al 



Kogut et al 



2005; Montrov et al 



2003)). Altogether, CMB tempera- 



ture and polarisation anisotropy measurements now enable us to constrain cosmological param- 
eters and hence the underlying theoretical models. A firm detection of the first peaks in the CMB 
anisotropy angular power spectrum and a damping at smaller scales has now been obtained. 
However, an excess of power at scales around a few arcminutes (I > 2000) has been found in the 
data (see Sect.Q. 

In the present study, we consider only a subset of data which allows us to explore the whole 
range of measured scales. Namely, we use WMAP and BOOMERanG data, whose major advan- 
tage is an accurate measurement of the primordial signal at t < 1500, and that of ACBAR, BIMA 
and CBI (at higher €) which exhibit an excess of power with respect to the expected primary con- 
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tribution. In addition, we take into account in our analysis the polarised data from WMAP and 
BOOMERanG. Altogether, this sample represents 1470 data points. 




10000 



Fig. 1. Data used in the analysis. Overplotted the best model when only primary fluctuations are 
considered and fitted to data with t < 1500 (purple solid line). 



3. Contributions to the signal 

Observations of the sky may contain many astrophysical processes which contribute to the signal 
along the line of sight. They have to be either accounted for or removed before analysing the 
signal of interest. Usually CMB experiments observe high latitude patches of the sky from which 
galactic contaminations are well subtracted. In the following, we will therefore neglect potential 
galactic residuals and focus exclusively on the extra-galactic contributions described below. 



3.1. Primary fluctuations 

We call primary fluctuations those anisotropies that are present at recombination. They may 
be due to intrinsic temperature fluctuations, density fluctuations or velocity fluctuations in the 
primordial plasma. The derived CMB angular power spectrum shows a particular shape the main 
characteristics of which are independent of the cosmology: A plateau at large scales, a succession 
of peaks and a damping at small scales. The details of those features are, on the contrary, highly 
dependent on the cosmological parameters. The position of the first acoustic peak is sensitive 
for example to the total density of the Universe. The CMB power spectrum is now well used to 
determine cosmological parameters by fitting the data of Fig. 1 with theoretical predictions. Here, 
the temperature and polarisation angu lar power sp ectra for primary fluctuations were computed 



using the Boltzmann code camb (ILewis et al. 2000) 
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3.2. Secondary fluctuations - SZ effect 

The Sunyaev-Zel'dovich (SZ) effect is actually the dominant secondary effect at scales between 
a few tens arcminutes to a fraction of arcminute. It has now been observed in the direction of 



a few tens of known clusters (e.g. 



Carlst rom et al.ll2.Q02l and references therein). The SZ effect 



has a particular spectral signature which induces an excess of brightness at high frequency and a 
deficit of brightness at low frequency in the Rayleigh- Jeans regime. This makes it in theory easily 
distinguishable from other contributions to the millimetre (mm) and sub-millimetre (submm) 
sky when multi-frequency observations are conducted. When the sky is observed at one single 
frequency (e.g. CBI, VSA, BIMA, AMIBA Lo et alJ 2001) it is a priori impossible to clean the 



signal from an SZ contamination except in the direction of known clusters. An SZ contribution 
is thus expected at small angular scales in CMB surveys. It is usually estimated through the SZ 
angular power spectrum which can be inf erred directly from cosmological numerical simulations 



or computed analytically (e.g. 



Komatsu & Seliak 2002 for a review, and references therein). 



In the present study, we choose to compute the SZ power spectrum using an analytic ap- 
proach. Such a choice is particularly important for a coherent analysis of the CMB anisotropies 
including primary and secondary signals. Indeed, the SZ power spectrum depends on the cos- 
mological parameters through the number of SZ sources and their characteristics (gas content, 
temperature, evolution...). Accounting for modifications through numerical simulations requires 
to run numerous simulations with different parameters which is both time and storage consum- 
ing. Numerical simulations are also limited by resolution and size. This affects the estimates of 
contributions at small and large angular scales. The dependence on cosmological paramet ers is 



more directly and easily handled with analytic computations (see for example 



Komatsu & Seliak 



20021) . In practise, we have computed the power spectrum of sources wit h masses between 



M ffiitt = 5xl0 u /r 1 andM, 



= 5 xl0 l5 h 1 solar masses up to z = 7 following 
(2002). The Poisson contribution to the SZ power spectrum C^ z is given by: 



Komatsu & Seliak 



C 



sz 



Jo dz JM m 



dN(M,z) 
dM 



[y e (M,z)fdM 



(1) 



where V(z) is the comoving volume element, dN(M,z)/dM is the comoving mass function and 
ye{M,z) is the 2D Fourier transform of the projected Compton parameter (the line of sight 
integral of the intra-cluster gas pressure) which characterises the amplitude of the SZ effect. 
The sp ectral dependence is given by g{v). We used a modified Press-Schechter mass function 



Sheth & Tormenll99 9s and accounted for possible variations of the spectral index /,\ of the ini 
tial density power spectrum through the "degenerated shape parameter" defined in IwillkoOll) . 



Komatsu & Seliak (2002) 



For the individual contribution ye{M, z), we followed the description in 
of a hot gas in hydrostatic equilibrium in a universal dark matter potential with a constant poly- 
tropic equation of state. Detail ed effects of the gas profile and mass function on the SZ power 



spectrum are given in 



Komatsu & Seliak 



(2002) to which we refer the reader. Using eq. we 
can compute the expected signal from a population of SZ sources for any given set of cosmolog- 
ical parameters, at every multipole I and for any observing frequency v. The SZ power spectrum 
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has a bell-like shape the amplitude of which depends mainly on the normalisation parameter 
cr 8 , but also on the mass and redshift ranges together with the set of physical and cosmological 
parameters used in the computation. A correlated contribution to eq. Q would modify the low 
I part of the spectrum by adding power at large angular scales. The relative contribution of SZ 
anisotropies at those scales being small we choose not to take the correlations into account. The 
EE and BB polarisation power spectra induced by clusters are orders of magnitude smaller than 
the primary polarised power spectra ( Liu et ali2Q05 ). We therefore do not take these into account 
in the present analysis. 



3.3. Point source contamination 



The extra-galactic contaminations depend on the observing frequency through the spectral en- 
ergy distribution of the sources and their number counts. The extra-galactic contributions to the 
CMB signal are usually classified broadly into radio galaxies emitting mainly at low frequencies 
and dusty IR galaxies at high freq uencies with a c rossing point at about 150 GHz where both 



contributions are of the same order ( Aghanim et al. 2005). 



Star forming galaxies emit, in the near-IR, optical and UV domains and also in the far-IR (up 
to 90% of their radiation) and submm, an emission associated with the radiation absorbed and re- 
emitted by dust. The total emission from unresolved and faint so-called IR galaxies is responsible 



2000 



Miville-Deschenes et al. 2002; 



for the Cosmic Infra- red Back ground, the fluctuat ions of which are now detec ted llLagache et al. 
iMatsuhara et alJkoOO 



Kiss et alJl2001l) . This is an impor- 



tant contaminant for the SZ and CMB measurements at high frequencies (> 100 GHz). 

The radio galaxy emission is dominated by synchrotron emission from relativistic elec- 
trons spiralling around magnetic field lines. The synchrotron spectrum generally follows a 
powe r law V, with -1 < a < -0.5. However, those spectra need not be simple power 



laws (Herbig & Readhead 



1992) but may have curved shapes instead; they even can exhibit 
inverted slopes (with a > 0); moreover some of the radio sources may be vari able. Radio 



sources are rath er well cat alogued at low frequencies thanks to NVSS (Condo n et all 19981) and 



FIRST (White 



( Benn ett et al 



et al. 



1997) surveys. The all sky survey from WMAP observations at 41 GHz 



2003) has provided us with information at higher frequencies where the radio 
sources are poorly known. In order to assess the level of contamination due to radio sources 
at the frequencies of interest for the CMB (between 30 and 200 G Hz), studies were con ducted 



and models propose d to predict radio source c ounts ( Toffolatti 



de Zotti et all2005l) . In particular, the model by 



de Zotti etal. (2005 



etal 



1998 



Sokasian et al. 2001 



) found to be in quite a good 



agreement with detections of radio sources by WMAP provides us with good estimates of the 
radio source contribution to CMB signal. 

Point source contamination can be reduced in, or removed from, CMB observations. This is 
usually done using multi-frequency data and component separation techniques, or using moni- 
tored observations at the same frequency as CMB experiments, or finally using spatial correla- 
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tions with already existing source catalogues. Different strategies for point source cleaning were 
used by the groups o f ACBAR, BIMA and CBI, an d are d escr ibed in detail in the arti cles by 



Dawson et al 



(2002); 



Mason et al 



(2003); 



Kuo et al. 



(2004and 



Readhead et al 



(2004a). These 



strategies comprise: Pointed observations at the CMB frequencies of sources at lower frequen- 
cies, direct counts using deep and mosaic images at the observing frequency, removal of faint 
sources from external catalogues, ancillary surveys at lower frequencies, using existing cata- 
logues. . . However, the uncertainties in the physical description of the contaminating emissions 
(extrapolations from frequency to frequency, source variability. . . ) together with the sensitivi- 
ties of the complementary data leave an unremoved contamination that is likely to contribute to 
the signal. Tof folatti et al.l ( 120051) have recently revisited the contributions of extra-galactic point 
sources to the arcminute scale measurements of ACBAR, BIMA and C BI at the relevant fre 
quencies, namely 28.5, 30 and 150 GH z. In particula r using predictions by 



updating the source model of 



ToffolattLsLal 



de Zotti et al. 



(2005) 



( 1998), they have shown that faint sources with 
fluxes too weak to be detected in ancillary surveys might contribute a significant fraction of the 
arcminute signal. 



For our study, we consider that such an unremoved component does indeed contribute to 
the measured arcminute signal. We do not take into account possible enhancement by lensing. 
We model the point source contribution in a very simple way, assuming it is well described 
by a power spectrum C^ s = a^^ 2 , typical of a Poisson distribution. The parameter o^ s ac- 
counts for the r.m.s. contribution of unremoved point sources (i.e. below the detection limit or 
unaccounted for) at a given observing frequency v b s . As a consequence, a different power spec- 
trum, i.e. a different a% s is assumed for each experiment (see Sect. 0J. Correlations between 
sources modify the contribution at large angular scales by adding power to the sho t noise and 
thus they are expected to c hange the slope of the associated power spectrum ((e.g. 



2003 



Negre llo et al 



Song et al 



2005). However, in view of the still limited knowledge we have of it, we 



choose not to take the source correlation into account. 



Finally, point sources may in principle also contribute to the polarised signal on small scales. 
Predic tions on C MB polarisation level from extra-galactic radio sources are summarised by 



Tucci et al 



( 2005), who show that radio galaxy contribution is largely dominated by synchrotron 
emission from our galaxy for all relevant frequencies. Similarly, dusty IR galaxies are also ex- 
pected to contribute to the polarised signal. Not much data is available on IR galaxy polarisation, 
and firm predictions on its contribution to the CMB are difficult to make. Assuming the level 
of dust pola rised emission in external galaxies is not higher than that of the Milky Way, it has 
been shown I Tucci et alJ2005 ) that polarisation from IR galaxies is sub-dominant at frequencies 
below 1000 GHz and for ( < 3000. We therefore neglect polarisation from radio and IR galaxies 
in our analysis. 
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4. Analysis 

4.1. Method 



In this work, we account for the contributions to CMB signal of the different astrophysical 
sources (primary fluctuations, secondary anisotropics and residual extra-galactic point sources) 
in a coherent way. In other words, we compute their associated power spectra by taking into 
account the dependences on the cosmological parameters of each signal. Note however that our 
empirical modelling of residual point sources does not explicitly exhibit such a dependence. We 
analyse the CMB data shown in Fig. 1 and look for the set of best cosmological parameters (best 
fit model) and corresponding error bars when assuming one, two or three of the contributions 
described in section[3]contribute to the measured signal. 

As seen in the previous section, the SZ angular power spectrum is frequency dependent. We 
thus consider three observing frequencies 150, 30 and 28.5 GHz. The level of contamination 
from residual point sources is related to the beam size of each experiment, to its sensitivity and 
to the flux limit which accounts also for possible source extraction using cross-correlations with 
other wavelengths. Therefore, the point source contribution in terms of power spectrum is fitted, 
in the following, independently for each experiment. As a result, the total power spectrum (sum 
of primary, secondary and unremoved point sources) is also frequency dependent and thus fitted 
for each observing frequency, i.e. for each small scale experiment (ACBAR, BIMA, CBI) added 
to WMAP and BOOMERanG. 

We use a Monte Carlo Markov Chain approach allowing us to account for a large number of 
parameters witho ut using int e nsive computational resources. We apply the test of convergence 



developed by 



Dunklev et al. 



(2005) to check the robustness of our results. The likelihood of 
WMAP data is computed by using the WMAP collaboration routines 1 whereas the likelihoods of 
the other data (AC BAR, BIMA, CBI) are computed with the Bond, Jaffe & Knox approximation 



( Bond et al 



2000) when the necessary information is provided (Gaussian-like approximation 
otherwise). When all the astrophysical contributions are taken into account the final likelihood £, 
writes: 

N 

£ = Y]-C i (C ( (i,v i ,0)\Cf s (i)) 

i 

where C° bs (i) is the angular power spectrum of an experiment i observing at frequency v,, 9 is 
the set of cosmological parameters and the total theoretical angular power spectrum is: 



CS, v„ 9) = C™(0) + Cf(v,, 9) + Cf (vO, 

with C£ s (v,-) = al^ s i 2 . Note that in the case of WMAP and BOOMERanG experiments, the 
C^ nmary (#) (and a fortiori the likelihood £j) includes polarisation signal. We consider throughout 
a flat cosmological model without contribution from neutrinos nor tensor modes. We assume a 



1 These routines are available at 
http://lambda.gsfc.nasa.gov/product/map/current/m_sw.cfm 
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power law initial power spectrum for the density perturbations with spectral index n s . The re- 
maining cosmological parameters are: The baryon budget, Q.bh 2 , the cold dark matter content, 
Q c( j m /! 2 , the Hubble constant, H , the normalisation factor <j% (given by the amplitude of fluc- 
tuation on a sphere of 8/i~' Mpc), and the Thomson scattering optical depth r (related to the 
reionisation history of the Universe), to which we add the amplitudes, a^ s , of the point source 
contributions for each experiment. 

4.2. Results 

In the following we consider first the case where primary and secondary anisotropies are present, 
then we add the contribution of unremoved point sources. We show in both cases the set of de- 
rived cosmological parameters and point out the effects of an additional, non primary contribution 
to the signal. 

4.2.1. CMB+SZ case 

Atmospheric, galactic and resolved point sources are usually well removed from CMB observa- 
tions; SZ effect has thus firstly been advocated to explain the excess of power detected at small 
scales. Previous analyses often decorrelate anisotropies at large scales, associated with primary 
signal, from those at small scales related to secondary signal. A primary power spectrum is thus 
generally fitted with data at I < 1500. It is added to a template SZ power spectrum (generally 
obtained from numerical simulations with fixed cosmological parameters). In such analyses, the 
amplitude of the SZ power spectrum is the only parameter set free (often labelled cr^ z and re- 
lated to the normalisation cr 8 ). Such a template spectrum is then fitted to the data available above 
( = 1500 and the SZ contribution, and therefore cr^ z , are derived. 

As seen in Sect. [3] the SZ angular power spectrum depends on the cosmological parameters 
through the volume element, the number density and the individual Compton parameter. The 
relative contribution of primary and secondary signals are thus strongly connected and cannot be 
formally separated. As a result, although the signal from primary fluctuations is weak at small 
scales (f > 1500) it is incorrect to decorrelate the two contributions in the fit. 

In the present study, we thus fit all data described in Sect.|2]by angular power spectra sum- 
ming the primary and secondary contributions, with free cosmological parameters. The results of 
the corresponding MCMC run (Run 1) are shown in terms of one dimensional likelihood func- 
tions (solid red lines) in Fig. [2] The latter are compared with the likelihood functions obtained 
by fitting Z < 1500 data by primary CMB power spectra only (black solid lines). 

The main effect of adding coherently the SZ signal is to allow the cosmological models to 
fit the small scale data. This implies that high values of cr 8 (median value: 0.97) and n s (median 
value: 1 .02) are preferred. Both parameters tend to amplify the power at small scale. In particular, 
as seen before, the amplitude of the SZ power spectrum depends on < tx. The higher <xs the 
larger the SZ signal. Such an effect was recognised in previous studies Jpoldstein et alJl2003t 



10 
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Fig. 2. Marginalised one-dimensional distributions for the cosmological parameters investigated. 
The black dotted-dashed line shows the curves for the case where data with £ < 1500 are fitted 
with pure CMB models. The red solid lines are derived from the MCMC run with all data and 
SZ contribution added to CMB (Run 1). The value of D./,h 2 as been multiplied by 100 for clarity. 



Bond el al]l2005l i. A coherent analysis taking into account both primary and SZ signals allow s 
to exhibit the effects of a small scale excess of power on the other cosmological parameters. 
As a matter of fact, in order to compensate the additional small scale power and keep the large 
scale power at WMAP amplitude, large values of the optical depth t (close to 0.3) are necessary 
to balance the effects of erg and n s . Such values are indeed allowed b y the WMAP TE power 
spectrum at large scales (see iKogut et al.ll2003t iMelchiorri et al I2OO5L This is the well known 
degeneracy between erg or n s and r illustrated in Fig. [3] The other cosmological parameters Qt,/i 2 , 
£l m h 2 and Hq are only very marginally affected by the excess of power at small scales, but our 
results seem to suggest that slightly higher values of Q m /i 2 and Ho are favoured. 

The best model, defined as the one with cosmological parameters being the median of 
the distributions (Fig. 0, is shown in Fig. [4] The corresponding cosmological parameters are: 
Q. b h 2 = 0.025, Q. cdm h 2 = 0.116, H = 75.33 km/s/Mpc, r = 0.22, n = 1.02, cr 8 = 0.97. 
This model slightly improves the goodness of fit as compared to the best pure primary spectrum 
without adding any other parameter. In Fig.0] the corresponding total angular power spectrum is 
plotted in black solid line whereas the primary and secondary contributions appear respectively in 
purple dotted-dashed line and orange long-dashed line. The difference in SZ amplitude between 
the two plots at 30 and 150 GHz is due to the frequency dependence giy) 2 in eq. Q. The relative 
contributions to the total power spectrum are shown as functions of ( in Fig- 13 the primary con- 
tribution is displayed in white and the secondary contribution is displayed in orange/grey. This 
figure gives the contribution of the SZ signal observed in microwave in a joint and coherent anal- 
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Fig. 3. Two-dimensional plot of the degeneracy between erg and t obtained from Run 1 . Coloured 
shaded regions mark the 1, 2, 3 <x contours from light to dark 



ysis where cosmological parameters are adjusted simultaneously for both primary and secondary 
anisotropics to fit current data. At t — 3000 the CMB signal is dominated by the SZ contribution 
which represents 86% and 59%, at 30 and 150 GHz respectively, of the total signal. As one can 
see in Fig[2] the best model that fits i < 1500 data, free from SZ contribution, is still a good model 
for the whole data set (up to £ = 10000). As a matter of fact, the SZ contribution at 30 (150) GHz 
at i — 3000 varies from to 98% (91%) for the models within the 68% confidence interval in 
the 6 dimensional cosmological parameter space. The actual data set at small scale is clearly 
poorly constraining due to the large error bars and the inter-calibration problem. The situation 
is then expected to improve with future observations. Furthermore, it is worth noting that the 
present analysis does not account for the statistical dispersion which induces a large dispersion 
in the SZ contributio n. Most of the SZ s ignal is indeed due to rare massive clusters. As noticed 
in Fig. 2 of < 



' Coorav & Melchiorri (2002), the SZ variance is expected to be much larger due to 
the non-Gaussian nature of the SZ signal than that of a Gaussian distribution. Mutli-frequency 
observations of the same field will then be the only way around this problem. 



4.2.2. CMB+SZ+PS case 

As seen in Fig. |4j the SZ angular power spectrum does not fully explain the excess power of 
some data points at small scales. Of course, the latter could be lcr events. Nevertheless, as dis- 
cussed above, some contribution from unresolved and unremoved point sources is also expected 
in the small scale signal. We want to quantify the effects of this additional contribution. We thus 
repeat our analysis by including the point sources in addition to primary and SZ contributions 
(MCMC Run 2). Modelling the angular power spectrum by 1(1 + 1)Q /2/tt = a$f (i/2000) 2 , we 
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Fig. 4. Total power spectrum (in solid black line) overplotted to the data for the two frequencies 
considered. The SZ power spectrum amplitudes are function of the frequency range considered 
(magenta dashed lines) while primary (purple dot-dashed line) contribution is not. 
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Fig. 5. Contributions of primary and SZ power spectra to the total spectrum (solid black line of 
Fig- El- The orange shaded (white) region gives the percentage of total power due to SZ (resp. 
primary fluctuations) contribution. 



introduce 3 new parameters aj, s noted A, B, C and describing the amplitudes of point sources re- 
maining in each experiment, ACBAR, BIMA and CBI respectively. These new parameters give 
the amplitude in /iK 2 of the point source contribution at t = 2000. 

The output of Run 2 is shown in Fig. [6] again in terms of one dimensional likelihood func- 
tions. We notice in this case that the estimated cosmological parameters (red/light lines) are all 
in perfect agreement with the results obtained by considering Z < 1500 data and primary fluctu- 
ations only (black/dark lines). The addition of an unremoved point source contribution at small 
scales improves the goodness of fit more than when adding SZ only. The values of cosmological 
parameters Q.hh 2 , Q m /i 2 and Ho are almost identical to those obtained with primary CMB only 
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and i < 1500 data. The parameters r, n s and <x 8 also converge to values obtained in pure primary 
CMB analyses. In the present case (primary CMB, SZ effect and unresolved point sources) we 
find a perfect agreement between large and small scale data converging towards the following 
median values for the cosmological parameters: Q/,/i 2 = 0.0237, Q. c d„,h 2 = 0.111, Ho = 74.6 
km/s/Mpc, t = 0.16, n = 0.99, erg = 0.87. 
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Fig. 6. One-dimensional distribution of cosmological and point sources amplitude parameters 
from MCMC Run 2. The vertical (blue) lines in the three lower panels give the amplitude of 
point source contributions as computed by Toffolatti et al. (see text). 



In addition to the cosmological parameters, we furthermore obtain estimations of the point 
source amplitudes for each experiment. The corresponding values are 181, 14, 170 fjK 2 for 
ACBAR, BIMA, CBI respecti vely. A recent study b ased on a physical model for point sources in 
radio and infra-red domains I Toffolatti et all2005 ) has shown that the excess of power at small 
angular scales could be explained alm ost entirely by the con tribution from unresolved point 



sources. We compare the predictions of 



Toffolatti et al 



( 2005) to our results given by the one 



dimensional likelihood functions of A, B and C in Fig. [6] We find that predictions from Toffolatti 
et al. are not only in agreement with our results but their values correspond to the maxima of the 
distributions in Fig.|6] 

We plot the contributions of each component to the signal in Fig.Qfor our best fit model. 
The relative contributions of the three astrophysical processes to the total signal are summarised 
in Fig. [H] The primary CMB is displayed in white, the SZ effect is in orange/light grey and the 
point sources are in red/dark grey. At ( — 3000, the relative contributions are 2, 40 and 9% for 
SZ contributions at (ACBAR, BIMA, CBI frequencies) and 91, 32 and 85% for the point source 
signals (at the same frequencies). Again, these numbers, corresponding to the best model, have to 
be taken with caution as they do not account for the larger statistical dispersions due to the non- 
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Fig. 7. Best model using all contributions. Primary in purple, SZ in orange, point sources in red. 




Fig. 8. Relative contributions of primary CMB, SZ effect and unremoved point sources power 
spectra. Primary CMB is displayed in white, SZ effect is in orange/light grey and point sources 
are in red/dark grey. From left to right, the vertical lines mark the {$ observed by CBI, BIMA, 
ACBAR experiments. 

Gaussian nature of both SZ and point sources. Furthermore when the 68% confidence interval in 
6 dimensions is considered, the SZ (PS) contributions at t = 3000 vary from to 80%, to 96% 
and to 91% (0 to 99%, to 85% and to 98%) for respectively ACBAR, BIMA and CBI. 



5. Discussion 



As we have seen in Sect. 14.2.11 the excess power on small scales could be attributed to the sole 
contribution of SZ sources (86% and 59% at 30 and 150 GHz at I = 3000). In that case, high 
values of <r 8 , n s and of the Thomson scattering optical depth t are favoured. This in turn may 
suggest that the Universe went through a very extended reion isation period, which re mains com- 



patible with present observational constraints as reported by 



Melchiorri et al 



(2005). However, 
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once we include the possibility of unremoved or unresolved point source contribution, the like- 
lihood functions of t and <x 8 that we obtain peak essentially at the same values as obtained with 
I < 1500 primary CMB data only. 

Our method parametrising, through the parameters A,B,C, the extra-galactic source power 
spectrum enables us to probe the possible contributions to the signal. Furthermore, it allows to 
test the consistency of specific PS physical models with the cosmological model. In case of a 
mismatch, our method would imply the need for different physical descriptions of the source 
(number counts, evolution, spectral energy distribution). 

A substantial fraction of point source contribution seems allowed (see Fig.[8]and Sect. 14.2.21 . 



Such a contribution is in agreement with extra-galactic number count predictions ( Toffolatti et al. 



2005). We noted however that statistical dispersion due to the non-Gaussian nature of SZ and 
point sources is likely to reduce significantly the actual contributions. In addition given the cur- 
rent data, large ranges of PS contribution are allowed. From the CMB standpoint, despite the 
excellent agreement of Toffolatti's model, on the one hand, the SZ effect could effectively con- 
tribute and solely explain the excess of power at small scales. On the other hand, different point 
source models are allowed by current data. 

Despite their quality, current data are still not constraining enough for such kind of multi- 
contribution studies. In addition to their intrinsic error bars, the use of multiple experiments 
introduces additional inter-calibration uncertainties. Ideally, we would need an experiment that 
would observe continuously from £ = 2 to ~ 10000. The Planck satellite will only observe up to 
£ ~ 2000. However, we will still lack data at high and intermediate multipoles, especially in the 
range £ e [2000, 10000]. In a few years from now, a combination of Planck and future ground 
based high resolution experiments will help us in reducing significantly the error bars and solving 
partly for the inter-calibration problem. 



6. Conclusions 

We presented the first coherent analysis of CMB data available from £ = 2 to i = 10000, 
taking simultaneously into account contributions from the primary fluctuations and secondary 
anisotropies as well as point source contamination. More specifically, we have explicitely taken 
into account in a consistent way the cosmological dependence of secondary anisotropies. In the 
first part of our work, this allowed us to confirm that SZ effect contribution to the CMB power 
spectrum could explain the so-called small scale power excess, as found in previous studies. We 
showed however that this requires not only high values of <x 8 , but also high values of the spectral 
index n s and of the Thomson scattering optical depth t, leaving other cosmological parameters 
basically unaffected. 

In the second part of our work, we included also unresolved/unremoved point source con- 
tamination. In that case, we recover the cosmological parameter values as determined from large 
scale CMB data only. Our consistent analysis introduces a parametrisation of the point source 
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signal. Consequently, we obtain estimates of the contribution of each CMB source to the total 
power spectrum. Our method hence allows us to test for cosmological consistency of physical 
models of point source populations. 
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